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Abstract. Chaotic dynamics with sensitive dependence on initial conditions 
Oh may result in exponential decay of correlation functions. We show that for 

one-dimensional interval maps the corresponding quantities, that is, Lyapunov 
exponents and exponential decay rates, are related. For piecewise linear ex- 
panding Markov maps observed via piecewise analytic functions we provide 
explicit bounds of the decay rate in terms of the Lyapunov exponent. In addi- 
tion, we comment on similar relations for general piecewise smooth expanding 
maps. 



1. Introduction 

Chaotic motion accounts on the one hand for the well-known phenomenon of 
sensitive dependence on initial conditions, that is, exponentially fast divergence of 
nearby orbits, and on the other hand for the phenomenon of decay of correlations 
or mixing. Both properties are intimately related with the observation that even 
low-dimensional chaotic systems share common features with random processes. 

This intuitive picture has been used as a basis to address some of the fundamental 
questions arising in nonequilibrium statistical mechanics [22, 23, 15]. In fact, it 
is a simple exercise to show that topological mixing implies sensitive dependence 
on initial conditions. At the measure-theoretical level, however, relating Lyapunov 
exponents, the quantitative measures for sensitive dependence on initial conditions, 
to decay rates of correlation functions is a more involved task. For instance, it is 
easy to construct simple maps with finite Lyapunov exponent and arbitrarily small 
correlation decay (see, for example, a Markov model described in [17]). 

Thus, at a quantitative level it is tempting to explore in some detail in which 
way the rate of correlation decay is linked with Lyapunov exponents, as both quan- 
tities are supposed to have a common origin. In more general terms and from a 
wider perspective this topic can be viewed as belonging to the realm of fluctuation 
dissipation relations for nonequilibrium dynamics, where one cause, the underly- 
ing detailed dynamical structure, is responsible for the approach to the stationary 
state, that is, for the decay of correlations, but, at the same time, is responsible for 
fluctuation properties at a microscopic level, or in our context for the sensitive de- 
pendence of initial conditions and a positive Lyapunov exponent. Furthermore, any 
relation between decay of correlations and Lyapunov exponents is of great practical 
interest, as the measurement of Lyapunov exponents, unlike the rate of correlation 
decay, is notoriously difficult to determine in real world experiments [1, 18]. In [31] 
it was even suggested to take correlation decay rates as a meaningful approximation 
for Lyapunov exponents. 

The problem we want to address can be illustrated by a very basic textbook 
example, probably considered for the first time more than two decades ago [2]. 

l 
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Consider a linear full branch map (see Figure 1) on the unit interval I = [0,1], 
that is, a map /:/—>/ having a finite partition of / into closed intervals Ik with 
pairwise disjoint interior such that (i) for each k, we have f(Ik) = I, and (ii) / 
has constant slope 7^ on each Ik- The physical invariant measure is given by the 
Lebesgue measure and the Lyapunov exponent with respect to this measure can be 
expressed in terms of the slopes 

A = £|I fc |ln| 7fc |, (1) 

k 

with \Ik\ = 1 / 1 Tfc I denoting the size of the interval Ik- 

The exponential rate of decay for correlation functions is determined by the 
negative logarithm of the second largest eigenvalue in modulus of the associated 
Perron- Frobenius operator. In this setting it is well known (see, for example, [27]) 
that eigenfunctions of the Perron-Frobenius operator are given by polynomials and 
that the corresponding eigenvalues v m can be expressed as 

v m = t^t— = 141— (m >o), 

with largest eigenvalue Xo = vq = 1 and the subleading eigenvalue Ai being the 
second largest in modulus, |Ai| = max{ | z/i | , 1*2 } . Thus, correlation functions of 
sufficiently smooth observables decay typically at an exponential rate a = — In |Ai|. 
Since |Ai| > ^2 > we obtain an upper bound for the decay rate a < — In which 
can now be related to the Lyapunov exponent (1). If we apply Jensen's inequality 
to the convex function ip(x) — — ln(x) we end up with 

a<-ln£|J fc |l<£|J fc |(-lnl)=2A. (2) 

k fk k Ik 

The estimate of the decay rate in (2) has been based on which contains positive 
terms only, even if the slopes have different signs. As a result the upper bound 
is given by twice the Lyapunov exponent. If all slopes have the same sign, say 
7fc > 1, then v x determines the subleading eigenvalue and the Lyapunov exponent 
itself yields an upper bound for the decay rate, that is, a < A. 




Figure 1. Diagrammatic view of a linear full branch map. 

In this article we will address the question to which extent this simple reasoning 
can be generalised to a larger class of systems. From a mathematical perspective 
there is a considerable amount of literature on the existence of invariant measures 
of one-dimensional maps, the corresponding Lyapunov exponent, estimates for the 
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spectra of the associated Perron-Frobenius operator, and a possible relation between 
Lyapunov exponents and the decay of correlation functions [12]. Expansiveness of 
the underlying map / : I — > I is to some extent the key ingredient to establish a 
physical measure /x, that is, an invariant measure which asymptotically characterises 
a large set of orbits, or in formal terms an invariant ergodic probability measure 
which is absolutely continuous with respect to Lebesque measure (see, for example, 
[25, 21, 3, 10]). 

The Lyapunov exponent for a Lebesque typical point x £ I is given by 

A = [ \n\f(x)\d^= lim -ln\(f n Y(x)\. (3) 

Regarding decay rates for correlation functions, a substantial part of the problem 
consists in giving the notion of the rate of correlation decay a proper meaning (see, 
for example, [12]). With the standard definition of a correlation function for two 
observables ip and ip, 

C v ^{n) = J^{f n (x))^(x)dfj,- J^(x)dfi J^{x)d(j,, (4) 

the exponential decay rate a v ^ governing the asymptotic behaviour of the corre- 
lation function may be formally introduced by 

a <p,ip — sup{s : lim sup \exp(sn)C v ^(n)\ < oo} . 

The rate a v ^ is sensitive to the choice of ip and ip and can be made arbitrarily 
small or large by special choices for the observables [13]. However, it is possible 
to define a rate of decay with respect to "typical" observables from some linear 
space H 

an = inf {a v ^ : ip, ip £ %} . 

We will call this quantity the mixing rate. Note that this rate still depends, for 
instance, on the degree of smoothness shared by the observables in %. 

One approach for determining bounds on the mixing rate an relies on refor- 
mulating (4) in terms of the Perron-Frobenius operator. The rate an is then 
determined by the subleading eigenvalue or by the essential spectral radius of the 
Perron-Frobenius operator. The desired relation between the Lyapunov exponent 
and the mixing rate expressed in (2) then becomes a lower bound for the subleading 
eigenvalue. There is a considerable body of literature on upper bounds for spectral 
values (see, for example, [25, 21, 10]) providing useful tools to establish ergodic 
properties of dynamical systems. However, to the best of our knowledge hardly any 
nontrivial lower bounds exist (see, however, [28], where an exponential lower bound 
for correlation functions of suspension semiflows is given). 

The main tool to establish the desired inequality consists in studying the prop- 
erties of the generalised Perron-Frobenius operator Cp with potential — /31n|/'|, 

^){x)= £ (5) 

For (3 — 1 this expression reduces to the Perron-Frobenius operator which has lead- 
ing eigenvalue one. The subleading eigenvalue Ai (or the essential spectral radius, 
if no subleading eigenvalue exists) determines the mixing rate an = — In | Ai | . In 
addition, the derivative of the largest eigenvalue vo((3) of (5) with respect to /3, 
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that is, the derivative of the topological pressure, determines the Lyapunov expo- 
nent (see, for example, [8] for a basic exposition). Finally, the required estimate 
follows from the convexity of the topological pressure and gives a lower bound for 
the subleading eigenvalue of the Perron-Frobenius operator. 

The main challenge is, of course, to put these ideas into practice. Restrict- 
ing to piecewise linear Markov maps considerably reduces the need to worry about 
mathematical subtleties, as the operator Cp admits finite-dimensional matrix repre- 
sentations when considering observables consisting of piecewise analytic functions. 
Thus, at a computational level all technical details reduce to straightforward matrix 
manipulations (see Section 2), allowing us to keep the presentation elementary and, 
at the same time, making the underlying ideas transparent. To keep the presenta- 
tion self-contained basic properties of piecewise linear Markov maps together with 
properties of Cp are summarised in the Appendix. In Section 3 we are addressing 
the question whether the bounds presented in Section 2 hold for general (non-linear) 
expanding Markov maps. We have compelling evidence that the estimate breaks 
down if analytic observables are considered (see also [11]), but the validity of our 
proposition can be restored if observables of bounded variation are considered. 



Let us consider a topologically mixing piecewise linear expanding Markov map 
/:/—>•/ with respect to the partition V = of the interval I (see 

Definition A.l for a formal account). Denote by fk := fi k the k-th branch of /. 
A nice property of such maps is that the corresponding transfer operators have 
finite matrix representations. In fact, the space of piecewise polynomial func- 
tions of degree less than M, that is, the space of functions of the form x i->- 
SfcLi Em=o a kmX m Xk(x), where \k is the characteristic function of the interval 
Ik, is an /^-invariant subspace. While it appears to be difficult to trace the earliest 
reference for this result, to the best of our knowledge the first experimentally rele- 
vant application appeared in the context of power spectra for intermittent dynamics 
(see [27, 30]). 

Now, in the natural basis of piecewise monomials, the operator Cp restricted to 
the invariant subspace mentioned above is represented by the (M + l) N x (M+1) N 
block upper triangular matrix 



of C\ on the space of piecewise constant functions (see, for example, [9, p. 176]), 
shows that the matrix elements of the block matrices T^ mn ^ {(3) are given in terms of 
the slopes 7^, the intercepts dk of the branches fk of the map /, and the topological 
transition matrix (Aki)i<k,i<N induced by / and V (see (16)) as follows: 



2. Piecewise linear Markov maps 
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The eigenvalues are determined by the diagonal blocks T^ mm \j3) with matrix ele- 
ments given by the first factor in (7). 

Note that since the underlying map is assumed to be topologically mixing, the 
matrices T^ mm \/3) are irreducible and aperiodic for even m. The Perron- Frobenius 
Theorem (see, for example, [16, p. 53] or [24, p.536]), now guarantees that T^ mm \j3) 
has a simple, positive eigenvalue, larger (in modulus) than all other eigenvalues, 
which we denote by v m {j3). Now, vo(/3) determines the topological pressure, given 
by P(/3) = ln^ (/3), which has the following well-known properties. 

Lemma 2.1. 

(i) P(1) = 0; 

(ii) P(f3) is a convex function of (3; 

(iii) P(f3) is analytic in (3; 

(iv) P'(l) = -A.' 

These properties have been established for a rather large class of dynamical 
systems (see, for example, [20]). In our context, however, they follow easily from 
the matrix representation above using elementary methods. As an example, we 
shall derive property (ii). To simplify notation, we shall denote the first block 
matrix T(°°\/3) by T(j3) for the remainder of the argument. Observing that P((3) = 
ln^o(/3) is obtained from InTr ((T((3)) n ) /n as n tends to infinity, convexity of the 
topological pressure, that is, 

P(i8i(l - t) + flat) < (1 - t)P{fix) + tP(fo) (0<t< 1), 

follows from a simple estimate using the Holder inequality 

ilnTr((T(A(l-t)+/? 2 t)n 



-In Yl T Wi(/Ml - *) + faWMil -t) + /3 2 t) ■ ■ ■ T ln _ lla (^(l -t)+ fat) 

lo,li,...,l n -i 



~ - 1" I E T loh(Pi) ■ ■■Ti n _ 1 i (Pi) I I T lah (fi 2 )---T ln _ llQ {p2) 

^lo,ll,...,ln — l J \lo-,h,---ln 

= (1 - t)- InTr ((T(/3i_)) n ) + t- InTr ((T(/3 2 ))") . 
n n 

Here we have used that (7) implies Tki(/3\+02) = 2Tw($i)^fcj($2) as well as Tki(t/3) = 
(Tki (/?))*• The other statements of the lemma can be derived in a similar way, or 
can be found in standard textbooks on ergodic theory. 

To establish a relation between the correlation decay, that is, between the eigen- 
values of the operator (5) for j3 — 1, and the Lyapunov exponent, note that the 
largest eigenvalue of the Perron- Frobenius operator is given by fo(l) = 1, while 
^2(1) is a positive eigenvalue, which provides a lower bound for the subleading 
eigenvalue of Cp. Thus, on the one hand 

an < -Ini^(l). (8) 

On the other hand T^(f3) = T^(/3 + 2) by (7), which implies 

v 2 (/3)=v (/3 + 2). (9) 
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Hence, using the properties of the topological pressure in Lemma 2.1, the relations 
(8) and (9) yield 

a H < -P(S) < (3 - 1)A . (10) 
See Figure 2 for a graphical illustration of this result. 



2A 
A 
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Figure 2. Schematic representation of the topological pressure 
P{(i) and graphical illustration of the estimates (10) and (11). 

Note that if all slopes 7^ of / have the same sign, then we have T( n )(/3) = 
sign(7fe)T( 00 ^(/3 + 1). Thus, we can apply the previous arguments to |fi(l)| to 
obtain the following improved estimate 

a n < -In |^(1)| =- In 14,(2) = -P(2) < (2-l)A. (11) 

To summarise, we have shown the following. 

Proposition 2.2. Let /:/—>/ denote a topologically mixing piecewise linear 
expanding Markov map. IfH — H(Dr) denotes a space of piecewise analytic ob- 
servables (see Appendix, Definition A. 2) then the mixing rate is bounded in terms of 
the Lyapunov exponent A with respect to the piecewise constant invariant density, 
by 

a n <2A. (12) 
If all slopes have the same sign the sharper estimate 

a n <A (13) 

holds. 

Remark 2.3. The assumption that / be topologically mixing is sufficient but not 
necessary. Indeed, there exist piecewise linear expanding Markov maps / with the 
following properties: the map / is not topologically mixing, yet exhibits exponential 
decay of correlations and the conclusions of Proposition 2.2 hold. 

The simplest examples for which the bounds are achieved are the tent map 
{an = 2A) and the doubling map {an = A). It is finally worth mentioning that 
the condition P{1) = 0, that is, the absolute continuity of the reference measure 
with respect to Lebesgue measure, is not essential for Proposition 2.2 to be valid. 
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The conclusions also hold, for instance, for certain Gibbs measures with respect to 
piecewise constant potentials. 

3. Remarks on general expanding Markov maps 

The setup of piecewise linear Markov maps is rather special. One may thus be 
tempted to ask whether a result like Proposition 2.2 extends, say, to Markov maps 
with finite curvature. The previous considerations are based on (9), a result which 
has been exploited previously in a slightly more restricted setup [2] to conjecture 
an exact relation between correlation decay and generalised Lyapunov exponents. 
It has already been demonstrated that such an identity breaks down for maps with 
finite curvature [11]. This fact, taken on its own, however, does not prevent the 
validity of a generalisation of Proposition 2.2, making a study of how curvature 
affects the previous considerations a worthy task. 

Using an Ulam-like construction any expanding Markov map can be approxi- 
mated by a piecewise linear map (see, for example, [4]). Let F : I —¥ I denote a 
piecewise expanding, but not necessarily linear Markov map with Markov partition 
{Ii, . . . , In}- Using the cylinder sets Ui 0i ... t i n _ 1 = {~^Z F~ k (Ii k ) we may introduce 
a piecewise linear approximation /„ : I — > I by the following construction. The map 

f n linearly interpolates F on each cylinder set, that is, ,f n {Ui a » n _i) = f/»i,...,» n _i 

(see Figure 3(a)). It is a straightforward exercise to show that the Lyapunov ex- 
ponent of /„ tends to the Lyapunov exponent of F (with respect to the absolutely 
continuous invariant measure) as n tends to infinity. However, the analysis of the 
mixing rate requires greater care. While Proposition 2.2 is still valid for any order 
n of the approximation, it is far from obvious whether the proposition is valid for 
a general expanding map. 

To illustrate this point we consider a simple example, a family of full branch 
piecewise Mobius maps F c defined on [—1,1], 

We restrict the parameter to c 6 (—1/4,1/2), in order to guarantee expansivity. 
Figure 3(a) depicts the map F c for c = —0.22. 

The leading part of the spectrum of the corresponding Perron- Frobenius operator 
considered on the space of analytic observables can be approximated using a spectral 
approximation method. The basic idea of this method is to approximate Cp by 
an n x n square matrix H n £pH n , where II „ denotes the projector that sends a 
function to its Lagrange-Chebyshev interpolating polynomial of degree n — 1. This 
method is easily implemented and, moreover, it is possible to show (see [5]) that 
the eigenvalues of H n CpH n converge exponentially fast to the eigenvalues of Cp. 
Using this method the leading eigenvalues of the Perron-Frobenius operator and 
their dependence on c are easily obtained (see Figure 3(b)). A minimum for the 
subleading eigenvalue occurs at about c = —0.11. The corresponding numerical 
value reads Ai w 0.10415 resulting in a mixing rate an = — ln|Ai| ~ 2.2619. 
The corresponding Lyapunov exponent (which hardly depends on the parameter 
c) is computed using the numerical approximation of the invariant density. The 
numerical value is A w 0.685 so that the inequality (12) is clearly violated. 

In order to understand why the reasoning at the beginning of this section fails, 
that is, why an approximation of the map F c by a piecewise linear Markov map 
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Figure 3. Left: (a) Diagrammatic view of the Mobius map (14) 
for c < (solid). Dotted lines indicate the cylinder sets of the 
first generations. Broken lines show the first two piecewise linear 
approximations /i and ji. 

Right: (b) The largest eigenvalues (in modulus) of the Perron- 
Frobenius operator for the map F c in dependence on c, as ob- 
tained using the Lagrange-Chebyshev method with truncation 
order n = 25. Positive/negative eigenvalues are indicated by 
filled/open symbols. 



fails to produce the correct mixing rate, let us consider increasingly finer piecewise 
linear approximations /„ of the map (14), see Figure 3(a). For each map /„ the 
Perron-Frobenius operator restricted to piecewise polynomial functions has a finite 
matrix representation (see (6)). For the remainder of the section we shall only deal 
with the case /? = 1 and refer to T((3) (and v m ({3)) as T (and v m ). Figure 4(a) 
shows the numerical results for the leading eigenvalues of the diagonal block T^ 11 ) 
and T^ 22 \ that is, v\ and v 2 , respectively, for increasing level of approximation 
n = 1, . . . , 6. For comparison we display the subleading eigenvalue v^ wb ^ of T^ 00 ' 
as well. For every level n of approximation, v 2 gives the subleading eigenvalue Ai 
of the Perron-Frobenius operator of the piecewise linear approximation, and these 
values seem to converge as n tends to infinity. The values and the limit are larger 
than exp(— 2A), meaning that the inequality (12) in the Proposition 2.2 is satisfied, 
as expected. The eigenfunction u n : [—1,1] — > R of the corresponding eigenvalue 
vi is a quadratic polynomial on each element of the partition, but it develops an 
increasing number of discontinuities between different intervals of the increasingly 
finer partition. These eigenfunctions do not seem to converge to a smooth limit 
(see Figure 4(b)). In fact, unlike the invariant density there is no reason why the 
limit should be smooth. The numerical experiment suggests that we end up with a 
function of bounded variation. 

It is indeed possible to show and perhaps well known that an estimate like (13) 
holds for discontinuous observables. For that purpose let us consider the Perron- 
Frobenius operator L\ on the space of functions of bounded variation. Recall that 
a function / : [— 1, 1] — \ K is of bounded variation if it has finite total variation 
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Figure 4. Left: (a) Leading eigenvalues v\ and V2 of the corre- 
sponding matrix blocks T^ 11 ) and T^ 22 ) for a piecewise linear ap- 
proximation of the map (14) with c = —0.11, as a function of the 
level of approximation. The subleading eigenvalue of T^ 00 ' , v^ uhl ^ , 
is displayed as well. For comparison, exp(— A) and exp(— 2A) are 
depicted as well. The broken lines are a guide for the eye. 
Right: (b) Eigenfunction u n corresponding to for n — 1, ... ,6 
with normalisation J I \u n (x)\dx = 1. For clarity, successive ap- 
proximations are shifted by 0.5. The open symbols indicate the 
discontinuity set of the eigenfunction, i.e., the increasingly finer 
Markov partition of the piecewise linear approximaton. 



var(/) = sup{£)i=i \f(xi) - /(aii-i)| : -1 < x < . . . < x p < 1} < oo. In this 
setup, the spectrum of the Perron-Frobenius operator associated with expanding 
maps has been studied in detail (see [19]). In particular, there is an explicit formula 
for the essential spectral radius given by 

a ess = lim (M{\(F*)'(x)\ : x e [-1, l]}T 1/k . 

A;— >oo 

Thus, we have an upper bound for the mixing rate 

a B v < -ln<T ess = lim - lninf{|(F c fe )'(a;) : x £ [-1,1]}, 

fe— >CX3 K 

which yields the following estimate for the Lyapunov exponent: 
k=\J\n\{F k c )\x)\dn 

> iinf{ln|(F c fc )'(x)| :xe [-1,1]} J d» 

= ilninf{|(^)'^)|:^6[-l,l]}. 
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Thus, for observables of bounded variation we have the following result. 

Proposition 3.1. Let f : I I be a piecewise monotonic smooth expanding inter- 
val map which is mixing with respect to its unique absolutely continuous invariant 
measure. Then the rate of decay of correlations for functions of bounded variation 
is bounded by the Lyapunov exponent 

a BV < A . (15) 

In fact, almost identical statements can be found in [12], for example, Corol- 
lary 9.2. 

4. Conclusion 

There is no simple, straightforward answer to the question about the relation 
between Lyapunov exponents and mixing rates. On formal grounds one may ar- 
gue that both quantities probe entirely different and independent aspects of a dy- 
namical system, and that no particular relation should be expected. Lyapunov 
exponents are determined by properties related to the largest eigenvalue of the 
Perron-Frobenius operator. By contrast, correlation decay depends crucially on 
properties of the observables, with mixing rates being related to the subleading 
part of the spectrum. Thus, abstract operator theory on its own does not seem 
to provide further insight into the relation between both quantities. Witness, for 
example, the doubling map viewed as an analytic map on the unit circle. While 
its Lyapunov exponent is finite, the Perron-Frobenius operator has no nontrivial 
eigenvalue when considered on the space of analytic functions, that is, correlations 
of analytic observables decay faster than any exponential (see, e.g., [3]). 

The argument outlined above, however, is a bit too simplistic. In fact, our re- 
sults on piecewise linear expanding Markov maps observed via piecewise analytic 
functions or general piecewise smooth expanding maps observed via functions of 
bounded variation suggest that bounds on the mixing rate in terms of Lyapunov 
exponents can be derived provided that specific properties of the underlying dy- 
namical system are taken into account. Estimates of this type rely on nontrivial 
lower bounds for spectra. As such, they are complementary to estimates which are 
available for proving the existence of spectral gaps, and will thus require completely 
different approaches. 

It turns out that the spirit of the result contained in Proposition 3.1 can be 
understood by considering an observable with a single discontinuity. In order to 
substantiate this claim we have performed numerical simulations on the map (14) for 
c = —0.11. We have computed the autocorrelation function (4) for the observable 
tp = ijj = 4>h with 4>h(x) = x if \x\ < 1/2 and 4>h{x) = x — siga(x)h if |a;| > 1/2, 
having a discontinuity of stepsize h at \x\ — 1/2 (see Figure 5). Choosing h = 0, the 
corresponding observable <ph is analytic and the correlation decay is seen to follow 
the subleading eigenvalue Ai of the Perron-Frobenius operator defined on the space 
of analytic functions (see Figure 3(b)). In the discontinuous case corresponding to 
h ^ the short time initial decay of the correlations still follows the pattern of the 
analytic observable, but the correlation function now develops an exponential tail 
which obeys (15). The tail becomes more pronounced if the stepsize increases. In 
fact, the mixing rate seems to be very close to the Lyapunov exponent. Revisiting 
the considerations leading to Proposition 3.1, it is tempting to surmise that this 
coincidence is a consequence of large deviation properties of finite time Lyapunov 
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exponents, since the expression for the essential spectral radius involves an extreme 
value of a finite time Lyapunov exponent. Thus, the relation between Lyapunov 
exponents and mixing rates for observables of bounded variation could be viewed 
to arise from the same mechanism already exploited in [2] for analysing the simple 
case mentioned in the introduction. 
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Figure 5. Normalised autocorrelation function of the discontinu- 
ous variable (f>h with different stepsizes h for the map F c (14) with 
c = —0.11. The straight lines represent an exponential decay with 
rate In |Ai| and — A. Ergodic averages have been computed as time 
averages of a series of length 2 x 10 4 for 5 x 10 8 uniformly dis- 
tributed initial conditions, skipping a transient of length 100. The 
horizontal dotted line indicates the order of magnitude of statisti- 
cal errors induced by the finite ensemble size. 

This simple demonstration gives support to the folklore that correlation decay 
is linked with Lyapunov exponents. Even if a real world phenomenon is sufficiently 
well-modelled by a smooth dynamical system for which to date no link between 
correlation decay and Lyapunov exponents can be established, one should keep 
in mind that modern data processing inevitably involves digital devices, which 
correspond to discontinuous observations. Therefore, in formal terms observables 
of bounded variation could be the relevant class for applications and in these cases 
Proposition 3.1 applies. 

At an intuitive level it is easy to understand why discontinuous observations 
result in correlation decay related to Lyapunov exponents. A discontinuous observ- 
able is able to distinguish between different "microstates" at a "macroscopic" level, 
that is, a discontinuous observation is able to distinguish two states at infinitesimal 
distance. As the distance between two nearby phase space points separated by 
a discontinuity grows according to the Lyapunov exponent, the sensitivity at the 
microscale may be transported to the macroscale, that is, it may filter through to 
the correlation function by a discontinuous observation. 
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In our context the mathematical challenge is to establish a relation between mix- 
ing rates and Lyapunov exponents for natural classes of observables, for example, 
full branch analytic interval maps observed via analytic functions. Besides the need 
for developing tools to obtain lower bounds for spectra, establishing the relation 
alluded to above also requires a deeper understanding of which dynamical feature 
causes the point spectrum of the Perron-Frobenius operator, that is, which signa- 
ture of the microscopic dynamics survives if viewed via analytic observables. This 
is reminiscent of coarse-graining approaches in statistical mechanics, for example, 
the introduction of collective coordinates and quasi-particles. Thus, tackling the 
mathematical problem above may well shed some light on some of the most fun- 
damental problems in contemporary nonequilibrium statistical physics of complex 
systems. 
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In this section we shall provide a short account of the technical details to make 
the results of Section 2 rigorous. The main thrust of the argument is to define 
a suitable function space on which the generalised Perron-Frobenius operator is 
compact. Results of this type for general analytic Markov maps are not new (see, for 
example, [29] or [26]). The special case of piecewise linear Markov maps discussed 
below, where a complete determination of the spectrum is possible, is probably 
well known to specialists in the field. Unfortunately, we are at loss to provide a 
reference for the results in Proposition A. 5 below, so we will outline a proof for the 
convenience of the reader. 

To set the scene we define what is meant by a piecewise linear Markov map. 
Before doing so we note that by a partition of a closed interval / we mean a finite 
collection of closed intervals {h,. ■ ■ ,In} with disjoint interiors, that is, int(ife) fl 
int(//) for k ^ I, such that UfeLi Ik = I- 

Definition A.l. An interval map / : I — > I is said to be a Markov map if there 
exists a finite partition {Ik}k=i °f I such that for any pair (k, I) either /(int(Jfc)) fl 
int(/;) = or int(/;) C /(int(Tfc)). If this is the case, the corresponding partition 
will be referred to as a Markov partition and the N x N matrix A given by 



will be called the topological transition matrix of the Markov map /. 

A Markov map / with Markov partition {Ik\k=i is said to be expanding if 
| f'(x) | > 1 for all x G int(Tfc). It is said to be piecewise linear if /' is constant on 
each element of the Markov partition, that is, f'(x) — 7& for all x G int(Tfc). 
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Appendix: Spectral properties of Markov maps 




otherwise 




(16) 



LYAPUNOV EXPONENTS AND EXPONENTIAL DECAY OF CORRELATIONS 



13 



Finally, we call an expanding Markov map with topological transition matrix 
A topologically mixing 1 if there is a positive integer p such that each entry of the 
matrix A p is strictly positive. 

In what follows we shall concentrate on topologically mixing piecewise linear 
expanding Markov maps. Our aim is to define suitable spaces of observables on 
which the associated generalised Perron- Frobenius operator or transfer operator 



(see (5)) is well defined and has nice spectral properties. It turns out that these 
spaces can be chosen from spaces of functions which are piecewise analytic. 
In order to define these spaces we require some more notation. 

Definition A. 2. Let D denote an open disk in the complex plane. 

(i) We write H°°(D) = { h : D -> C : h holomorphic and sup zeD \h(z)\ < oo } 
to denote the space of bounded holomorphic functions on D. This is a 
Banach space when equipped with the norm \\h\\ jjooim = sw P z ed IM^OI- 

(ii) We use H(D) = H ca (D) to denote the space of TV-tuples (hi,..., Hn) 
of bounded holomorphic functions on D. This is a Banach space when 
equipped with the norm ll/iH^^) = max{||/ifc|| ff oo(£)) : k — 1, . . . ,N}. 

The desired space of observables will now be defined by linking the disk D 
occurring in the definition above to the dynamical system as follows. Given a 
piecewise linear expanding Markov map with Markov partition {Ik}k=i> ^ us 
denote by tpki ■ Ik — > h the inverse branch of the Markov map from partition 
element Ik into the partition element 7; as well as its obvious analytic continuation 
to the complex plane. Observe now that, since the map is expanding, all inverse 
branches are contractions. We can thus choose two concentric disks D r and Dr of 
radius r > and R> r, respectively such that 



It turns out that H(Dr) is a suitable space of observables for the map, in the 
sense that the associated transfer operator (17) is a well defined bounded operator 
on H(Dr). This is the content of the following result. 

Proposition A. 3. Given a piecewise linear expanding Markov map f with topo- 
logical transition matrix A and inverse branches ipki, suppose that the disk Dr is 
chosen as above. Then, for any real (3, the transfer operator Lp is a well defined 
bounded operator from H(Dr) into itself and is given by 



Proof. The representation (19) follows from a short calculation using the definition 
(17) of Cp. Since ^^(z)^ = |7;|~^ is constant and the disk Dr satisfies (18), the 



This is a slight abuse of terminology, since its use is usually restricted to continuous maps. 
However, it serves the same purpose as in the continuous setup as it guarantees the existence of 
a spectral gap for the corresponding transfer operator (see Corollary A. 6). 




(17) 



<Pm(Dr) C D r for all inverse branches ipkl ■ 



(18) 




(19) 
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operator maps T-L{Dr) to H(Dr). In order to see that Cp : W.(Dr) — > H(Dr) is 
bounded observe that if h G V.(Dr) with ||/i||^(£> H ) < 1, then 

W C P h W-H(D R ) = m ^\\(^ h )k\\ H oo {DR ) < max^] A ;fe | 7; |-^ < oo. □ 

Remark A. 4. The space %{Dr) is not the only suitable space of observables. Re- 
stricting to one and the same disk of analyticity for each branch, however, simplifies 
notation. More general spaces are discussed in [6] and [7]. 

Next we shall explain why the spectrum of Cp viewed as an operator on TI(Dr) 
is given by the eigenvalues of the diagonal blocks T( mm '(/3) for m G No defined in 
(6). 

The key ingredient of the proof of this statement is a factorisation of the transfer 
operator together with an approximation argument. 

In order to explain the factorisation of the transfer operator we observe that in 
(19) the argument of hi, that is Lpki{z), is contained in the smaller disk D r because of 
(18). We can thus use (19) to define the operator on a larger function space, namely 
H(D r ). Note that the space is 'larger' as analyticity is only required on a smaller 
disk D r C Dr. We shall write Cp : H(D r ) — > %{Dr) in order to distinguish this 
lifted operator from the one occurring in Proposition A. 3. Note that the arguments 
in this proposition can be adapted to show that Cp is a bounded operator. 

It is tempting to think of Cp and Cp as being essentially the same, since they 
are given by the same functional expression. However, as operators the two are 
different as the latter is defined on a larger domain. Yet, both operators are related 
by restriction. In order to give a precise formulation of this fact we introduce a 
bounded embedding operator J which maps the smaller space %(D R ) injectively 
into the larger space %(D r ). To be precise J : %(Dr) — > H(D r ) is given by 
(Jh) k = Jh k , where J : H°°(D R ) -> H°°(D r ) in turn is given by {Jh)(z) = h(z) 
for z e D r . Note that J looks superficially like the identity. This, however, is 
misleading as argument and image are considered in different spaces. 

The relation between Cp and Cp can now be written as 

Cp = CpJ. (20) 

Note that the factorisation above disentangles the intricacies of the map contained 
in Lp from its general expansiveness contained in J . 

We now turn to the approximation argument. For piecewise linear Markov maps 
the transfer operator is easily seen to map piecewise polynomial functions of degree 
at most M into piecewise polynomial functions of degree at most M. This follows 
from a straightforward calculation using the fact that the inverse branches are affinc 
functions. 

In order to exploit this property of the transfer operator further we shall in- 
troduce a projection operator defined as follows: given an analytic function h in 
H°°(Dr) and an integer M we use Puh to denote the truncated Taylor series 
expansion 

(p M h)(z)=f: f ^(z-z ) k , 

where z denotes the centre of the disk Dr. Clearly, Pm is a projection operator. 
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It turns out that the projections Pm approximate the embedding J for large M 
in a strong sense. In order to make this statement, the heart of the approximation 
argument alluded to above, more precise, we observe that, by Cauchy's Integral 
Theorem, we have for any h £ H°°(D R ) and any z G D r 

«*> -<**><«>- si i^l^f^*- 

where the contour Y is the positively oriented boundary of a disk centred at zq with 
radius lying strictly between r and R. It follows that the norm of J — JPm viewed 
as an operator from H°°(D R ) to H°°(D r ) satisfies 

R / r \M+l 

\\J - JPm\\ h ~ {Dr) ^ H co (Dt) < j^— r (^J 

In particular, we have 

JjmJJ- JPM\\ H oo (DR) ^ H oo (Dr) =0. (21) 

In order to extend this result to the space of piecewise analytic functions H(D R ) 
we introduce the projection operator Vm ■ T~L{Dr) — > 'H{Dr) by setting Vuh — 
(Puh\, . . . , Puhjsi)- The analogue of (21) now reads 

Jm^ \\J JV M \\ H{DR) ^ n{Dr) = . (22) 

We are now able to combine the factorisation (20) with the approximation result 
above to prove the main result of this Appendix. 

Proposition A. 5. Suppose we are given a piecewise linear expanding Markov map 
f with inverse branches ipkl an d disks D r C Dr satisfying (18). Then, for any real 
/3, the transfer operator Cp viewed as an operator on H(-Dfl) is compact and its 
non-zero eigenvalues (with multiplicities) are precisely the non-zero eigenvalues of 
the matrices T^™" 1 ^/?) with m £ N given in (7). 

Proof. We start by recalling that for every M > the transfer operator Cp leaves 
the space V M {U{D R )) invariant, that is, Cp{P M {'H{D R ))) C V M (H(D R )). Thus 

(I-V M )C fi V M = 0, 

where I denotes the identity on %{D R ). 

Using the above equation and the factorisation (20) we see that 

\\Cfj -VmCpPm\\ U ( Dr ) = Pm)\\ h{Dr) 

= \\CfsJ(l -V M )\\n(D R ) < \\£fi\\n(D r )^n(D n )\\J - JPhi\\n(D R )^n(D r ) , 
which, using (22), implies 

a)Too ~ V m^m\\ u{Dr) = . (23) 

Since PmZp'Pm is a finite-rank operator for every M, the limit above implies that 
Cp is compact. Clearly, the non-zero eigenvalues of each Vm^-^Pm are exactly the 
non-zero eigenvalues of the block matrices (6). The remaining assertion, namely 
that the non-zero spectrum of the transfer operator is captured by the non-zero 
spectra of the finite dimensional matrix representations follows from (23) together 
with an abstract spectral approximation result (see [14, XI.9.5]). □ 

Specialising to topologically mixing Markov maps we obtain the following refine- 
ment of the above proposition. 
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Corollary A. 6. Suppose that the hypotheses of the previous proposition hold. If 
the Markov map f is also topologically mixing, then Cp : %{Dr) — > T-L(Dr) has a 
simple positive leading eigenvalue Moreover, this leading eigenvalue is the 

Perron eigenvalue of the matrix T^ 00 ) (/?) . 

Proof. This follows from the previous proposition together with the observation 
that for m > 1 the spectral radius r{T^ mm \j3)) of the matrix T^ mm \p) is strictly 
smaller than the Perron eigenvalue of T^ m \j3). In order to see this note that for 
all m > 1 we have 



\Tt m \P)\<CT^\f3), 



where 



C = 



< 1. 



inf; 1 7; 

A short calculation shows that for each k > and each m > 1 we have 



where ||-|| F denotes the Frobenius norm. The spectral radius formula now implies 
that 

r{T^ mm \l3)) < Cv {(3). D 
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